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We present a new signature by which to one could potentially discriminate between a spectrum 
of gravitational radiation generated by a self-ordering scalar field vs that of inflation, specifically 
a comparison of the magnitude of a flat spectrum at frequencies probed by future direct detection 
experiments to the magnitude of a possible polarization signal in the Cosmic Microwave Background 
(CMB) radiation. In the process we clarify several issues related to the proper calculation of such 
modes, focusing on the effect of post-horizon-crossing evolution. 



It was pointed out over 15 years ago that a continu- 
ous symmetry breaking phase transition, can produce a 
scale invariant spectrum of gravitational radiation sim- 
ilar to that of Inflation pQ. However, until prospects 
arose, via measurements of polarization of the CMB, 
for an actual measurement of this background, no de- 
tailed quantitative estimate of this effect had been per- 
formed. More recently a quantitative calculation was 
carried out in the case of a scalar field that breaks a 
global O(N) symmetry [2|, confirming qualitative fea- 
tures and also demonstrating that the phase transition 
signal could be significantly larger than previously ex- 
pected . 

Since the detection of primordial gravitational radia- 
tion is widely considered a smoking gun for inflation it 
is important to carefully explore how to separate and 
disentangle this potentially competitive signal. Bau- 
mann and Zaldarriaga [5] have recently confirmed that 
that the problem is more subtle than one might ini- 
tially suppose. The standard representation of CMB 
polarization data in terms of spin-weighted harmonics 
of the B-modc would not be able to discriminate gravi- 
tational radiation produced by Inflation from that pro- 
duced by a phase transition. Instead, in order to probe 
for superhorizon-sized correlations, they introduced a 
new kind of i?-mode dubbed the B mode in ref [3] and 
have argued that angular correlations on scales slightly 
larger than 2 degrees could distinguish an inflationary 
signature from that due to other phase transitions . 
Subsequently Adshead and Lim [Jl [5] have pointed out 
that non-Gaussianity may provide another diagnostic. 
Both of these signatures will be challenging to extract 
from CMB data, however. 

The purpose of this short note is is to demonstrate an 
in-principle clear and direct distinction between gravi- 
tational waves generated as a result of Inflation vs self- 
ordering, which derives from the fact that a scalar field 
associated with symmetry breaking continues to act as 
a source of gravitational waves inside the horizon in the 



latter case. Primordial gravitational waves produce a 
CMB polarization signal due to anisotropic scattering 
of radiation by electrons in a quadrupole field, and this 
effect is dominated by modes which enter the horizon 
at decoupling, and again during reionization, both well 
into the matter-dominated era. However, all currently 
proposed and envisaged direct interferomctric gravita- 
tional wave detectors such as LISA and its projected 
successors will be sensitive to modes that entered the 
horizon in the radiation dominated era. As we shall 
demonstrate, there will be differing relative signatures 
for CMB signals vs those that might be observed in 
LISA or its successors in each of the two different cases 
(inflationary produced waves and those resulting from 
self-ordering scalar fields). 

In order to calculate this difference one must care- 
fully follow the dynamics of the scalar field inside the 
horizon. In so doing we demonstrate the need to calcu- 
late effects well inside horizon crossing . We also clarify 
confusion in the literature regarding the normalization 
conventions appropriate when comparing a CMB po- 
larization signal and a direct sub-horizon gravitational 
strain signal, and discuss the effects of vacuum energy 
domination on the evolution of the signals. 

Like all forms of radiation, gravitational waves be- 
gin to oscillate as their period approaches the age of 
the universe and they come inside the horizon. At the 
same time, they begin to redshift. For gravitational 
waves from inflation, which are only sourced during 
the inflationary epoch, these are the only effects that 
determine how the initial frequency-independent ini- 
tial strain value translates into a gravitational wave 
amplitude on horizon scales. This is reflected in the 
behavior of the strain power spectrum as a function of 
momentum wave vector; in the inflationary case dur- 
ing the matter dominated era the spectrum goes as 
3ji(fcr)/(fcr) whereas in the radiation dominated era 
the dependence goes as jo, (where jo and j\ are spher- 
ical Bessel functions, k is the comoving wavenumber, 
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and r is the conformal time). These are equal mag- 
nitude as fcr — > 0, but differ by 15% at horizon cross- 
ing (with a larger amplitude during matter domination 
compared to radiation domination) (i.e. see [6]). 

In order to properly calculate the magnitude of 
waves inside the horizon one must take into account 
not just redshifting but also the possibility of addi- 
tional sourcing. In the case of self-ordering, waves are 
continually sourced on the horizon scale and below be- 
cause fields begin ordering on this scale. Indeed, this is 
what is responsible for the generation of horizon-sized 
modes in the first place [U [2] . This process continues 
for some period after horizon crossing, and given the 
differing dynamics during these periods there is no rea- 
son to expect that the results will be the same during 
these two epochs. 

In the existing literature a number of approximations 
are made for ease of calculation; these must be handled 
carefully in order to properly estimate this effect of the 
sourced radiation, and it is to this task we now turn. 
In [2] the evolution of the scale factor was smoothly 
interpolated between radiation dominated and matter 
dominated eras, and the evolution of modes were fol- 
lowed from their production (which was assumed to 
occur in the matter dominated era- relevant to CMB 
polarization) to well inside the horizon. In later work, 
Fenu et al [7] focused on waves entering the horizon 
during radiation domination. 

To investigate the effect of post-horizon crossing 
growth we plot the evolution of such a mode from its 
production at time tq to a time when its frequency is an 
order of magnitude smaller than that of the horizon, 
calculated by numerically solving the coupled scalar 
field evolution equations along with the metric during 
radiation domination. As can be seen in Figure 1, not 
taking this subsequent evolution into account reduces 
the results by almost an order of magnitude. (The 
oscillatory behavior presumably reflects damped ring- 
ing of the wave modes upon sudden termination of the 
source generating them.) 

Using these detailed numerical estimates for the 
strain power generated during these two epochs we 
can compare the contribution to the energy density to- 
day in modes which crossed the horizon during radia- 
tion domination vs matter domination for backgrounds 
from self-ordering scalar fields and from Inflation, when 
these are normalized to give identical contributions in 
either the high frequency or low frequency domain to- 
day. We assume, for the purposes of simple compari- 
son, an exactly flat inflationary spectrum here. 

As can be seen in Figure 2, sourcing of modes in- 
side the horizon is clearly more effective during mat- 
ter domination than radiation domination, effectively 
combatting redshift damping for some time. This may 
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FIG. 1: Appropriately normalized power in Fourier modes 
with p = xr as a function of pr. The solid curve represents 
a calculation including full scalar field dynamics down to 
pr — 10, while the dotted curve cuts of sourcing at horizon 
crossing, pr = 1. 

be due to the fact that during matter domination the 
energy momentum tensor of the relaxing scalar field 
differs more dramatically from that of the background 
matter density, therefore more effectively coupling the 
changing stress tensor to the production of gravita- 
tional waves. In any case if direct detection can be ac- 
complished with a precision of ±50%, then the distinc- 
tion between waves generated by cosmic phase tran- 
sitions and those by inflation should be easily distin- 
guishable. In practice of course, the spectral tilt of the 
gravitational wave spectrum would have to be mea- 
sured in the CMB in order to extrapolate expectations 
to frequencies appropriate for gravitational wave in- 
terferometers. In this regard it is also important to 
note that if waves from cosmic phase transitions were 
incorrectly attributed to inflation one would infer an 
incorrect frequency dependent tilt dependence for the 
assumed inflationary modes. 

Inflationary plots such as figure 2 are generated by 
numerical solution of the source- free gravitational wave 
equation. Weinberg has derived a number of asymp- 
totic analytic solutions to this equation that are useful 
for checking numerical work and for semi-analytic esti- 
mates [HI . The formulae as derived apply in the eras of 
radiation and matter domination. Since the universe 
is now entering an epoch of vacuum domination it is 
worthwhile to consider the extension of Weinberg's for- 
mulae past matter-vacuum equivalence. In the eras of 
interest the relation between conformal time and scale 
factor is given by the Friedman equation 

t = ^-( da' {n R + n M a' + n A a' 4 y 1/2 (1) 

where H n is the present value of the Hubble constant, 
and fli is the density parameter of the i'th form of 
energy. It is convenient to define a cq = (Qr/Qm) 
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component corresponds to modes that are pushed out- 
side the horizon and freeze; D\, to modes that continue 
to decay forever. Eq ^ shows that the eventual fate 
of the modes is a sensitive oscillatory function of the 
wave- vector q with period goo. 

Similar formulae can be written for modes that enter 
during the era of matter domination (q\ q e q)- 
The WKB solution in this case is 



FIG. 2: The fractional energy density per logarithmic fre- 
quency interval as a function of frequency (normalized to 
the frequency of modes that cross the horizon at matter- 
radiation equality) for gravitational wave modes generated 
during inflation compared to that in modes generated by 
scalar field reordering, for transitions tuned to give the same 
gravitational wave amplitude for high frequency modes that 
enter the horizon during radiation domination. 



and a,A = (Om /^a) 1 ^ 3 , the scale factors at radiation- 
matter and matter-vacuum equivalence respectively. 
It is also helpful to define the wave-vectors q cq — 
a eq H cq /\/2, q A = H A a\/V2 and q^ = l/r^. It 
is easy to see that the Friedman eq ([TJ) simplifies to 
t = 2[y/l + y — l]/<7 cq ,with y = a/a eq , in the eras of 
radiation and matter domination (a <C a\) and to 
r = — l/^A-z hi the era of vacuum domination 
a > a A . 

Weinberg showed that modes that enter the horizon 
in the era of radiation domination (q 3> q C q) are given 
by the WKB approximate formula 



Dg(r) = D q (T )— sinqr. 

ny 
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where k = q/q cq and Z3 g (ro) is the amplitude of 
the mode when it is far outside the horizon. The 
WKB formula applies once the mode enters the hori- 
zon (y ^> and we find it applies not only in the era 
of radiation and matter domination but it also contin- 
ues to apply through the present epoch deep into the 
era of vacuum domination until redshift z <C Q where 
Q = q/q\- We can continue the WKB solution eq ^ 
into the indefinite future, z 3> Q, by connecting it to 
exact solutions to the gravitational wave equation that 
apply in the era of vacuum domination, a 3> a\. We 
find 
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This formula applies for all z ^> 1. Here D\(x) — 
sin a; — a; cos a; and D 2 (x) = — xsinx — cos a;. The D 2 
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D q (r) = -D q (T )—^- cos qr 



(4) 



and gets extended to 
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Finally, to compare CMB limits and direct detection 
experiments, one needs to take into account different 
normalization conventions. CMB limits usually quote 
constraints on the strain directly, whereas traditionally 
constraints on direct signatures employ the energy den- 
sity per logarithmic frequency interval. A consistent 
convention is to write the metric as — a 2 (7]^ +h^) 
where 77^" is the Minkowski metric, a is the scale factor, 
and h denotes the strain perturbation. In this conven- 
tion we define the Fourier transform of the strain by 



where 



£ ^/<K00* A (k,T)(r) t 
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and the strain power in the modes is defined via 
<5 3 (k-k')P(fc,r)= £ <h* + (k,r)h + (k',r) 

+ h* x (k,T)h x (k',T) >,(8) 

where "+" and "x" label the polarization states, and 
the gravitational energy density normalized to the crit- 
ical energy density in the universe is defined via 
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din A: w gw (fc, r n ) 



(9) 



where oj gw (fc) is the normalized energy density per log 
frequency interval and r n is the present conformal age 
of the Universe. From the textbook formula relating 
the strain amplitude to the energy density of a gravi- 
tational plane wave it follows 
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is the relationship between strain power and energy 
density. 

An alternative convention [3] is to define the Fourier 
transform with a factor of (27r) 3 / 2 in the measure. This 
convention produces strain Fourier amplitudes that are 
smaller by a factor of (2-7r) 3 / 2 and a strain power that 
is smaller by a factor of (2n) 3 . The energy density 
per log frequency interval is still related to the strain 



power by eq ( 10 ) and it is therefore also smaller than 



the corresponding quantity in the first convention by 
a factor of (27r) 3 . As a result, the relation between 
the total gravitational energy density and the energy 
density per log frequency interval is then eq ([9| with 
the factor of (27r) 3 omitted. 

Ref [2] followed the first convention but with the 
twist that h + and h x differ from the above definition by 
a factor of —2. As a result the strain power and energy 
density calculated in that paper should be multiplied 
by a factor of 4 to be consistent with the convention 
we define here. In that paper, however a comparison 
of the gravitational wave energy density produced in 
this case with that due to inflation at horizon cross- 
ing used a quoted result for the latter based not on 
the same convention used to derive results for cosmic 
phase transitions but rather that of [9] . As a result the 
statement there that the ratio of the energy densities 
in the two cases is larger than the naive dimensional 
estimate (rj/vN 1 / 4 ) 4 by more than four orders of mag- 
nitude (where rj is the scale of the cosmic phase tran- 
sition and v 4 is the energy density during inflation) 
should be modified by a factor of 4/(27r) 3 = (27r) -3 
(i.e. see [7]). The correct statement is that the ra- 
tio of the gravitational energy density due to cosmic 
phase transitions to the inflationary energy density at 



frequencies accessible to direct detection is larger than 
the naive dimensional estimate by over two orders of 
magnitude. 



The results here, which now allow a consistent and 
detailed comparison of strain power and associated en- 
ergy densities a stochastic background of gravitational 
waves produced by cosmic phase transitions and infla- 
tion in two different regimes-for waves which entered 
the horizon during matter vs radiation domination- 
demonstrate the potential importance of comparing 
the magnitudes of the strain power in these two 
regimes. As we have shown, if such a background is 
observed via a CMB polarization signal, a comparison 
with a possible direct gravitational wave detection at 
shorter wavelengths would provide a powerful tool in 
order to distinguish the origin of such a background, 
and thus provide an important diagnostic to investi- 
gate possible new GUT scale physics. 
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